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^ ! Abstract 

■ We present a numerical study of the fermion-induced effective action in the pres- 

ence of a static inhomogeneous magnetic field for both 3 + 1 and 2 + 1 dimensional 
C<*) \ QED using a novel approach. This approach is appropriate for cylindrically symmetric 

magnetic fields with finite magnetic flux <3?. We consider families of magnetic fields, 
dependent on two parameters, a typical value B m for the field and a typical range d. 
We investigate the behavior of the effective action for three distinct cases: 1) keeping 
<J? (or B m d 2 ) constant and varying d, 2) keeping B m constant and varying d and 3) 
keeping d constant and varying $ (or B m d 2 ). We note an interesting difference in 
the limit d — > +00 (case 2) between smooth and discontinuous magnetic fields. In 
the strong field limit (case 3) we also derive an explicit asymptotic formula for the 
3 + 1 dimensional action. We study the stability of the magnetic field and we show 
that magnetic fields of the type we examine remain unstable, even in the presence of 
the fermions. In the appropriate regions we check our numerical results against the 
Schwinger formula (constant magnetic field), the derivative expansion and the numer- 
ical work of M. Bordag and K. Kirsten. The role of the Landau levels for the effective 
action, and the appearance of metastable states for large magnetic flux, are discussed 
in an appendix. 



1 Introduction 

The evaluation of the fermion-induced effective action in the presence of a static magnetic 
field is important in quantum electrodynamics. Historically, the effective action for a homo- 
geneous magnetic field was first computed by Heisenberg and Euler 0, then by Weisskopf 
0, and later by Schwinger |J with the proper time method. The effective action in this case 
can be expressed by an explicit formula. In reference Q Redlich obtained the corresponding 
results for the 2 + 1 dimensional QED. 

For inhomogeneous magnetic fields explicit formulas are not available. Approximations 
can be developed for the case of a) weak magnetic fields (perturbative expansion) and b) 
smooth magnetic fields (derivative expansion). 
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More recently, for a special form of the magnetic field it has been possible to derive a 
closed formula for the effective action for both 3 + 1 and 2 + 1 dimensions || ^j. Another 
exact result is presented in reference 0, for a magnetic field with a delta function profile in 
two dimensional Euclidean space. 

The aim of this paper is to obtain, by numerical computations, qualitative knowledge of 
the dependence of the fermion-induced effective action on various quantities characterizing 
the magnetic field such as its size and its degree of inhomogenity etc. For this numerical 
study we use a simple novel approach for the fast computation of the effective action, which 
is appropriate for static, cylindrically symmetric magnetic fields with finite magnetic flux $. 
We will call such fields magnetic flux tubes. 

Previous works towards the numerical study of the effective action are the following: 
in reference || the effective energy, in 3 + 1 dimensions, is computed numerically for a 
cylindrically symmetric static magnetic field which is constant inside a cylinder and zero 
outside it. In a more recent work || the effective action was computed for a cylindrically 
symmetric static magnetic field with a delta function profile, in 3+1 dimensions. Note that 
this field has infinite classical energy. 

Our purpose in this paper is to go beyond the previous numerical works, and to study, 
for the first time, the 2+1 QED effective energy in the presence of a magnetic flux tube. In 
addition, in the case of 3 + 1 dimensions we consider more realistic magnetic field configura- 
tions than those of the previous works, for example the Gaussian magnetic field of Eq. (2) 
bellow. 

We study families of magnetic fields dependent on two parameters, a typical value B m [] 
for the field (or the magnetic flux $ = B m d 2 /2) and a typical range d. We investigate the 
behavior of the effective action for three distinct cases: 1) keeping $ (or B m d 2 ) constant and 
varying d, 2) keeping B m constant and varying d and 3) keeping d constant and varying $ 
(or B m d 2 ). The mass of the fermion my is assumed to be constant. 

We note a difference in the limit d — > +oo (case 2) between smooth and discontinuous 
magnetic fields (section 5). We will see that this difference is related to the fact that the 
derivative expansion fails for discontinuous magnetic fields like the magnetic field of Eq. (3) 
below. 

An interesting problem not covered by previous works, is the behavior of the effective 
energy for large magnetic flux </> (strong magnetic field) and fixed rrifd, where <fi = e$/27r. In 
section 7 we present an investigation for large <fi. Especially in the case of 3 + 1 dimensions, 
we derive an explicit asymptotic formula for ^> 1. 

We also consider the question of stability of the magnetic field in the presence of the 
fermions. The question is, whether some radius exists that minimizes the total energy 
(classical energy plus effective energy) for fixed magnetic flux $, thus rendering the magnetic 
field stable at the quantum level. The same question was considered by the authors of 
reference [[| , for the magnetic field of Eq. (3) (below) in the case of 3 + 1 dimensions, and 
the answer was negative. Our purpose here is to study the stability of magnetic flux tubes 
in the case of 2 + 1 dimensions. In addition in the case of 3 + 1 dimensions, we generalize 
the results of reference p] for more realistic magnetic field configurations like the Gaussian 
magnetic field of Eq. (3). Note, that we answer this question for several two-parameter 

1 We define this typical magnetic field strength as B m = J B ■ dS/ird 2 . 
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families of magnetic field configurations of the type of flux tube for both 2+1 and 3+1 
dimensions (see the magnetic fields at the end of the present section). 

In appendix C we explain the role of the Landau levels for the effective action, relating 
them to the metastable states which appear in the case of large magnetic flux <3> (see appendix 
B). 

The vector potential and the magnetic field for the class of magnetic fields we examine, 
are given by the following equations: 

Ar) = \ f(r) e% B{r) = l|-(r 2 f(r)) e z (1) 

where /(r) is a function with asymptotic behavior 20/r 2 as r tends to infinity and = e$/27r. 
In what follows, we assume the rescaling B — > eB except where we state otherwise. We 
present our numerical results for two magnetic fields, one with a Gaussian profile and one 
which is constant inside and vanishes abruptly outside a cylinder of radius d with respective 
field strengths 

Bl(r) = | exp(-l!) (2) 
B 2 (r) = ^0(d-r) (3) 

where 6{x) is the step function. 

We have performed analogous calculations for six two-parameter families of magnetic 
fields other than those of Eqs. (2) and (3), with magnetic field strengths 

B s(r) = ^ —^-—^ , S 4 (r) = -| - exp(--), B 5 (r) = 



d 2 cosh 2 (r 2 /rf 2 ) ' 4W d 2 \d 2 ) FV d^' 5W rf 2 vr(r 4 /rf 4 + l) 
B 6 (r) = ^(l- r -)6(d-r), B 7 (r) = ^ (1 - ^) 6{d - r) , B 8 (r) 2 ' 



d 2 v d' v " ,w d 2 v rf 2; v " ow rf 2 (r 2 /rf 2 + l) 5 

Since including the results in detail would make this paper unnecessarily long it is perhaps 
sufficient to state that our conclusions (about the stability of the magnetic field or the 
dependence on d of the effective action etc) for the magnetic fields of Eqs. (2) and (3) are 
valid also for the above-mentioned magnetic fields. 

2 Effective energy 

We consider the following path integral that corresponds to QED 

Z = J DADTOfe'/^^-i^^+^W-™/)*) (4) 

in dimension D (D = 3,4). Integrating out the fermionic degrees of freedom in the above 
mentioned path integral, we obtain the effective action expressed as the logarithm of a 
determinant 

S eff [A] = - In f vyVVe i f dDx *W- m f )iB = -tr\n(tfl-m f ) (5) 
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where Ift = ^(d^ — ieA^), and 7 = 03, 7 1 = iai, 7 2 = 102 is a two component representation 
of the gamma matrices for the 2 + 1 dimensional QED. A four component representation 
will be used for the 3+1 dimensional case. 
We will use the identity 

tr\n(i fl - m f ) = ^tr \n(J7> 2 + m 2 f ) (6) 

in order to take advantage of the diagonal form of the operator lj) 2 +m 2 . It is well known that 
for 2 + 1 dimensional QED a parity violating term (Chern-Simons term) is induced [|], |i~0 



However, this term is zero for the case we investigate f\. See also the relevant comments in 



reference 1 1 



We deal first with the 2 + 1 dimensional problem. The operator p + mj for a static 
magnetic field can be put into the form 

2 + m 2 f = d 2 o + ( 1 rn D m ) 2 + m 2 f (m = l,2) (7) 

This operator has a complete system of eigenfunctions of the form tyi n y(x, t) = e~ lut ^ { n }{x) 
where ^{ n }(x) satisfies the eigenvalue equation (■j m D m ) 2 ^/{ n j(x) = i?{„}\I/{ ra }(x) and {n} is a 
set of quantum numbers. Note that E{ n } > since the operator (7 m Z) m ) 2 is positive definite 
as the square of a hermitian operator. It is convenient here to refer to the effective energy 
E e ff = —S e ff/T instead of the effective action S e ff, where T is the total length of time. 
The effective energy, if we perform a Wick rotation of the integration variable uj (uj — > iu), 
is given by 

I r+00 

#e//(2+i) = -— / du\n(uj 2 + E {n} + m 2 ) (8) 
Ztx {n} Jo 

We make the above integral convergent by subtracting a term which is independent of 
the magnetic field. We may do so because our purpose is to compute the difference between 
the effective energy in the presence of the magnetic field, and the energy when the magnetic 
field is absent. If we subtract the term — ^ S{n} Jo + °° dojhiiui 2 + m 2 ) from (8) and integrate, 
we get 

Eeff '(2+1) = -- (\/ E M + m / ~ m f) ( 9 ) 

M 

This expression incorporates the expectation that very massive fermions interact weakly with 
the magnetic flux tube (i.e the effective energy tends to zero as mf tends to infinity). 

The extension to the 3+1 dimensional case can be made if we replace m 2 by k 2 + m 2 in 
(9) and integrate over the k 3 momentum. An overall factor of 2 (from the Dirac trace) must 
also be included. For details see e.g reference 0. 

IT + E ^ + m ) ~ \Ai + m}) (10) 

-00 ZvT N v / 

where L z is the length of the space box in the z direction. 



2 The induced Chern-Simons term is of the form -| J d 3 x e^ vp A^ d v A p and it is zero, for our case, since 
A a = and A x = A 2 = 0. 
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Introducing a momentum cut-off A/2 we find: 



Hi) 

where we have omitted L z and we assume that we calculate energy per unit length. 

For large to/ the above expression for the effective energy must tend to the first diagram 
(vacuum polarization diagram) of the perturbative expansion of the effective energy: 

1 f +oc \~, m2 .m 2 f + q 2 x(l — x) 



£e//(3 + i) = ^3 J dxx(l-x) J o dqq\B(q)\ 2 \n(^ 



m) 



A 2 



+ < 12 > 

where B(q) = J d 2 x e~ tq ' x B(x). Note that B(q) depends only on q = \q\ as we have assumed 
that the magnetic field is cylindrically symmetric. 

If we compare the equations (11) and (12) for To/ — > oo we find 

Y,E {n} = -±-Jd 2 xB 2 (x) (13) 
M 

Thus, the logarithmically divergent part in equation (11) can be dropped if we include an 
appropriate counterterm in the QED Lagrangian. This counterterm is defined by the on-shell 
renormalization condition f\ 11(0) = 0. The renormalized effective energy is then given by 
the equation 

{n} / {n} 

Note, that the renormalized effective energy vanishes for large values of to/, as expected. 
For the sake of simplicity we will drop the index (ren) in the rest of this paper. 

In order to define the effective action for massless fermions for which the right hand side 
of Eq. (14) is singular, we should impose the renormalization condition for the vacuum 
polarization funtion Il(g 2 ) at a spacelike momentum — M 2 (IT(— M 2 ) = 0). Accordingly, a 
formula for the effective action for massless fermions is given in appendix A. 

In cylindrical coordinates and for the special magnetic field of Eq. (1), the operator 
{^j m D m ) 2 has the following diagonal form 

( 7 m D m ) 2 = -Dl-D 2 2 -sB 

Q2 i q i q2 i q r 2 f 2 (r) 

-f(r)+ V -sB(r) (15) 



For the definition of the vacuum polarization function n(<7 2 ) and for details about the renormalization 
condition for II (q 2 ) see, for example, the reference [12|. 
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where s takes the values ±1 that correspond to the two possible spin states of the electron. 
Because the operator {^( m D m ) 2 is invariant under rotations around the z axis, the eigenfunc- 
tions ^ have the form $?i n y(x) = e d< ^u(r)/ \/r, where I is the quantum number of the 
angular momentum (I = 0, ±1, ±2, ...). The function u(r) satisfies the equation 




k 2 u{r) (16) 



where 



v ljS (r) = -lf(r)-sB(r) + r -^l (17) 

and k 2 is the corresponding eigenvalue (we may set E{ n y = k 2 since E{ n y > 0). The spectrum 
of the radial Schrodinger equation (16) besides the continuous spectrum (0, +oo), also may 
have zero modes according to the Aharonov-Casher theorem 0|TB[. Note that for our special 
case {n} = {k, I, s}. 



3 Density of states 

In order to sum up over the continuous modes of the equation (16), we shall need the density 
of states % = pi, s (k) (i.e the number of states per unit k). 

puk)=tf: ee \k)+\ d -^ (is) 

where p/ s (k) = f is the density of states for free space, and Si yS (k) is the phase shift which 
corresponds to I th partial wave with momentum k and spin s. The relation (18) is shown e.g 
in reference |L4| . 



A noteworthy feature of the function Si }(T (k) is that it exhibits jumps which corresponds 
to met ast able states for large values of 0. This is shown in Fig. [| of appendix A. 

If we use the relations (18) and (9), the effective energy for 2 + 1 dimensional QED is 

£ e //(2+i) = ~ [°° (v^ 2 + ™? - m,) ^ (j: S l)S (k) j dk (19) 

where p/ s (k) was dropped since it contributes a field independent term. Note that the 
zero modes do not contribute explicitly in (19). 

The series of the phase shifts over I in (19) is not absolutely convergent. The simplest 
way to define it, is to sum symmetrically over I. We define the function 

L 

A(k)= lim Yl M*0 (20) 

L ^ +0 °s,l=-L 

The asymptotic behavior of the phase shifts for \l\ 3> kd is 

M*)-(-^ + fV (21) 



'According to this theorem a zero mode exists if the conditions 1 < I + 1 < <fi and s — 1 are satisfied. 
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and the sum 5i jS (k) + 5-i jS (h) vanishes for enough large values of I so the limit in (20) exists. 

The large I asymptotic behavior can be obtained from the WKB approximation for the 
phase shifts: 

M*) = ife | / + °° (jv-vUr)-^ ~ ^ 2 "^) ^ } ( 22 ) 

For large / only the asymptotic 0(1/ r 2 ) tail of the potential contributes to the above integral. 
Replacing v^ s (r) by its asymptotic form ^ ~ 2 2< ^ yields (21). 

Note that an asymmetric sum would give the same result since \im L _ t+OQ Yl^t=-L &i,s{k) = 
A(k) + d and d is independent of k. 

In appendix B we describe two ways for the numerical calculation of the phase shifts 
by solving an ordinary differential equation. The numerical calculation shows that A (A;) 
for large k tends to a constant value c = — 7r0 2 . An analytical proof for this result can be 
obtained using the WKB approximation (22) for the phase shift. In Fig. [I] we have plotted 
the function A(k) — c for = 4.5 and d — 1. The numerical computation was performed 
for the magnetic fields B\ and B 2 of Eqs. (2) and (3) respectively. In order to achieve 
convergence of (20) for = 4.5 we summed up to l max = 20 for the magnetic field B 2 and to 
Imax = 30 for the Gaussian magnetic field B\. 

Fig. |l] shows that the function A(k) — c tends rapidly to zero, apparently like l/k A for 
the magnetic field B 2 and faster for the magnetic field Bi, so the integral (19) is highly 
convergent. This is remarkable because it means that the same function A(k) can be used 
for the 3 + 1 dimensional case, as can be seen from the equation (26). 

The result that the function A(k) — c which corresponds to the magnetic field B\ tends 
faster to zero than the magnetic field B 2 suggests that in general the function A(k) — c 
spreads when the corresponding magnetic field is getting more localized. This has been 
checked by calculations for magnetic fields other than B\ and B 2 (see section 1). In addition 
we observed that when the magnetic field is getting more localized we need less values of I 
in order to achieve convergence of (20). 

The value of the function A(k) at k — is independent on the magnetic field configuration 
and depends only on <p. This becomes evident if we take into account the Levinson theorem 
which is presented in appendix A. Thus the two curves in Fig. ^ intersect at k — 0. 

Now from the relations (19) and (20) we obtain 

E eff{2+1) = -I l +C ° (v/P+^J - m f ) dk (23) 

Integrating by parts we find 

1 r+00 u 

Eefm + i) = 2- / , a (A(fc) - c)dk (24) 

An J o Jk 2 + mj 

For the 3 + 1 dimensional case, from the relations (14), (18) and (20), again dropping the 
free part from the density of states, we find 

1 /' +00 /, 2 2m f k2 + m 2 f \ d(A(k) - c) „ 

1 r-w)-^ (25) 



47r 2 Jo dk 
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Figure 1: A (k) — c as a function of k for <fi = 4.5 and d — 1 



Integrating by parts we obtain 



e//(3+l) 



2tt 2 Jo 



/ j^ 2 _i_ m 2 \ 

fcln I ;r-L {A(k)-c)dk 



m 



f 



(26) 



4 Dependence on the range d for fixed magnetic flux 

A dimensional analysis shows that A(k) is a function of the form F(kd, <fi). If we make the 
change of variable y = kd and set A(k) = F(kd, <p), we find 



E 



e//(2+l) 



e//(3+l) 



27rd Jo 
1 



{F{y,<f>)-c) dy 



y 2 + m 2 d 2 

+oo f y 2 + m 2 f d 2 \ 

V In ^ m2 J — j [F{y, (j)) - c) dy 



(27) 
(28) 



where d is a typical range of the magnetic field (see Eqs. (2) and (3)). 

In Fig. we see that E e ff( 2 +i), for the Gaussian magnetic field B%, is a positive decreasing 
function of d for fixed <fi. The asymptotic behavior of E e ff( 2 +i) for small values of d is 
proportional to 1/d. This can be seen from relation (27). For large values of d we are in 
the weak field regime and E e ffM+i) is given approximately from the vacuum polarization 
diagram 



E 



(2) _ 

e//(2+l) - g7r 



\ pi r+oo 

— / dxxil — x) I dy 
r z d Jo Jo 



y\B(y/d)\'< 



y (rrifd) 2 + x(l — x)y 2 



(29) 
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Figure 2: (ft 2 dE e ff( 2 +i) as a function of d for 
Gaussian magnetic field of Eq. (2). 



2.5,4.5,6.5,8.5 and mj = 1 for the 



where we have made the change of variable q = y/d and B(q) is the Fourier transform of 
the magnetic field. From the relation (29) we find that E e ffM+i) is proportional to 1/d 2 for 
large values of d. 

The curves in Fig. ^| appear to approach each other as d tends to infinity. Thus the 
effective energy for large d is proportional to 2 . This is expected because the effective 
energy for large d is given approximately by the vacuum polarization diagram which is 
proportional to <fi 2 . 

In Fig. |3| we present our numerical results for the magnetic field B 2 . We see that the 
effective energy as a function of d for fixed has the same features as those of the Gaussian 
magnetic field of Eq. (2) [|. Also we can see that the effective energy of B 2 is always bigger 
than the effective energy of B\ for the same <j) and d (this is also true for the corresponding 
classical energies). This fact seems to be a consequence of the flux of the magnetic field B 2 
being more concentrated than the flux of B\. 

The 3 + 1 dimensional case, for the case of B 2 , has already been studied with a different 
method by M. Bordag and K. Kirsten in reference ||. Our results are shown in Fig. f|. They 
agree closely with those of Fig. 3 of reference || . 

The 3 + 1 dimensional effective energy, for large d, is given by 



E 



(2) 

c//(3+l) 



8n 3 d 2 



dxx(l 



x 



dyy\B(y/d)\ 2 ln( 



m 2 d 2 + y 2 x(l — x) 
m 2 d 2 



(30) 



5 We have checked that this is also true for the 3+1 dimensional case. So the effective energy, in 3 + 1 
dimensions, is a negative increasing function of d for fixed magnetic flux (as shown in Fig. |j) for all the 
magnetic fields of the form of the magnetic field of Eq. (1). 
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Figure 3: (p 2 dE e ff( 2+ \) as a function of d for 
magnetic field of Eq. (3). 



2.5,4.5,6.5,8.5 and rrif = 1 for the 



For d — > +oo we find 



^e//(3+l) 



240vr 3 m2rf 4 j 



<% 3 |%/d)f 



(31) 



For further discussion we will need the Fourier transforms of the magnetic fields of Eqs 
(2) and (3): 

q 2 d\ 



Bi(q) 
B 2 {q) 



27r0exp( 



4 



47T0 



qd 



(32) 
(33) 



From the relation (31) we see that for large d the effective energy is proportional to l/d A . 
This is true only for magnetic fields for which the integral (31) over y is convergent. An 
example is the magnetic field B\ with the Gaussian profile. For the magnetic field B 2 the 
integral (31) is divergent. We found numerically from (30) that in this case the asymptotic 
behavior is 0(l/d a ), where 3 < a < 4. 

For small values of d, from the relations (28) and (13) we find the following asymptotic 
formula 



E, 



e//(3+l) 



12n 2 rrifd' 

The asymptotic formula (34) is the same as the one derived in reference 
(34) we obtain E eff{3+1) = ^ for tne magnetic field B 2 and # e //( 3 +i) 

for the magnetic field By. 



(34) 



From the relation 



6nd 2 
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Figure 4: 2 d 2 E e ff( 3+ i) as a function of d for = 0.4, 5, 10 and rrif = 1 for the magnetic 
field of Eq. (3). 

The total energy for B 2 is given by 

Etot(3+l) — -^dass(3+l) + ^e// (3+1) 

= 5^ + ^ /(m ^ 0) (35) 

where we have set £ , e //( 3+ i) = jpf(mfd, <fi). The physical meaning of the total energy is the 
energy we should spend in order to create the magnetic flux tube. In this article we have 
computed this energy at one loop order (i.e we have not taken into account virtual photons). 
This approximation is valid only for small values of the coupling constant e. 
For small values of d we find 



2vt0 2 2 1 
e 2 d 2 3nd 2 rrifd' 



Etot(3+i) = ~ ttt^ ln(— -) (36) 



We see that the logarithm, which comes from the effective energy part, dominates for small 
values of d. This asymptotic behavior of the total energy is not reliable. The reason is that 
small values of d correspond to magnetic fields whose Fourier transforms have contributions 
from high momenta at which the running coupling constant can not be assumed small. 
The asymptotic relation (36) for the total energy can be written in the form 

2tt0 2 

Etot{3+1) ~ ^W 2 (37) 
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where we have defined a d-dependent effective coupling constant, for small values of d, 
according to 



1 



1 



-e// 



(d) 



Gtt 



1 M 1 



rrifd 



(38) 



It is interesting to compare this to the B-dependent effective coupling constant of reference 
|I5|, for a homogeneous magnetic field, which for large B reads 



ln ( — ) 

/ 



12tt 



m 



(39) 



5 Dependence on the range d for fixed magnetic field 
strength B m = 20/ d 2 



0.4 



0.2 



<- 



'our numerical results' o 
'0.829' -— 
0.829-0.1 33/d 




Figure 5: E e ff( 2 +i)/d 2 as a function of d for B m = 5 and rrij = 1 for the magnetic field of 
Eq. (3). The discrete points correspond to our numerical results and the continuous line to 
the best fit curve of the form jd. The asymptotic value ao = 0.829 ± 0.002 agrees 

quite well with that calculated from the Schwinger formula: 0.829. Our numerical results 
for a and a\ are rounded to three decimal points. 

It is convenient to define a characteristic magnetic field strength as B m = J B ■ dS/nd 2 , 
in order to estimate the intensity of the magnetic field, where d is the range of the magnetic 
field. In the previous section we examined the dependence on d of the effective energy for 
fixed magnetic flux <fi. Here, we study the dependence on d (or irtfd) for fixed magnetic field 
strength B m = 2<p/d 2 (or B m /m 2 ). 
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Figure 6: E e ff^ + i)/d 2 as a function of d for B m = 5 and mj = 1 for the magnetic field of 
Eq. (3). The discrete points correspond to our numerical results and the continuous line to 
the best fit curve of the form jd. The asymptotic value ao = 0.129 ± 0.002 agrees 

quite well with that calculated from the Schwinger formula: 0.130. Our numerical results 
for a and eti are rounded to three decimal points. 



First, we examine the magnetic field of Eq. (3) (the Gaussian magnetic field of Eq. (2) is 
examined in the next section). An interesting feature of this field (of Eq. (3)), is that when 
its range d increases it tends to a homogeneous magnetic field. Thus, in the homogeneous 
limit (d ^> 1 / y/B m ) , we expect our numerical results to approach the result calculated from 
the Schwinger formula , which for the cases of 2 + 1 and 3 + 1 dimensions reads 



rpSch 

tj eff(2+l) 



A 



£3/2 f+00 fa 



in 



3/2 



,3/2 



coth(s) 



1 x > -sm 2 f /B 



- I e 

s 



(40) 



rpSch 

■^eff {3+1) 



A f 2 T ^( coth(s )_i_£)e--?/* 
871^ Jo s V s 3/ 



(41) 



where A is the area of the space box on the x — y plane and B is the magnetic field strength 
of the homogeneous magnetic field. 

In Fig. [| we have plotted the effective energy divided by d 2 for B m = 5 and rrif = 1 for 
the magnetic field of Eq. (3). We see that E e ff( 2 +i)/d 2 is a positive increasing function of d 
which tends very slowly to an asymptotic value. The determination of this asymptotic value 
is not possible with a directly numerical computation, as for large values of d the numerical 
error of our method becomes significant. However, we see in Fig. [5] that our numerical results 
(for d > 0.4) lie on a curve of the form ao — al/d. Fitting a curve of this form to our data we 
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obtain a = 0.829 ± 0.002 and a x = 0.133 ± 0.001. The asymptotic value a = 0.829 ± 0.002 
agrees quite well with that calculated from the Schwinger formula: E^j^ 2+1 ^/d 2 = 0.829. 
This asymptotic behavior has been also observed for other values of B m and rrif which cover 
all the range of the characteristic ratio B m /m 2 , that determines how strong is the magnetic 
field. 

For small values of d, if we take into account the vacuum polarization diagram of Eq. (29) 
which is good approximation for the effective energy, we find that E e ff (2+1) / d 2 is proportional 
to d. This is in agreement with the linear behavior of E e ff( 2 +i)/d 2 , for d — > 0, which is seen 
in Figs. and [7]. 

In the case of 3 + 1 dimensions, we see in Fig. |B| that our numerical results are negative 
and tend slowly to an asymptotic value, which also agrees closely with that calculated from 
the Schwinger formula of Eq. (41). The asymptotic behavior for large d is the same as that 
observed in the case of 2 + 1 dimensions (do — ai/d). 

For small d we obtain from the vacuum polarization diagram of Eq. (30) that E e ff (2+1) / d 2 
is proportional to In d. We see that this logarithmic behavior is not clear in Fig. [7] since we 
have not plotted enough points for small values of d. We avoided this, as in the region of 
small d (d ~ 0.1) the numerical accuracy we could achieve was not satisfactory. 

A different way, of finding the behavior of the effective energy for small d, is to replace 
F(y,cf)) = <pF^ l \y) + cf) 2 F^ 2 \y) in Eqs. (27) and (28) and take into account that the linear 
term is identically zero (F^(y) = 0). We can prove that = usine; the first Born 

approximation for the phase shifts. 

Our numerical work shows there is a significant possibility that the 1/d term in the 
above-mentioned asymptotic behavior ( QiQ — d\ /d) is true. If, indeed, this behavior is true, 
it would be interesting if one could gave an analytical proof. 

Note, that this asymptotic behavior is expected only for discontinuous magnetic fields 
like the one of Eq. (3). In Fig. we see that in the case of the Gaussian magnetic field of 
Eq. (2), our results for large d have the same asymptotic behavior with that of the derivative 
expansion: b§ — h\jd 2 . This means that there is a remarkable difference between smooth and 
discontinuous magnetic fields, in the way their effective energies (divided by d 2 ) tend to their 
asymptotic values. It is obvious that this difference, is related to the fact that the derivative 
expansion fails for discontinuous magnetic fields, like the magnetic field of Eq. (3). 



6 Derivative expansion 

An approximate way to compute the effective action in the presence of a static magnetic 
field, is the derivative expansion. This method should give accurate results for slowly varying 
magnetic fields. Our aim is to examine cases in which the derivative expansion fails to 
provide a good approximation. Also, we can test our numerical results comparing with 
those of derivative expansion in the homogeneous limit (B m d 2 ^> 1), where the derivative 
expansion is expected to give accurate results. 

The derivative expansion for 3 + 1 and 2 + 1 dimensional QED has been investigated in 
references [ 17|, [16| respectively. For the case of 2 + 1 dimensions, if we keep the first two terms 



(a zeroth order term plus a first order correction) the effective energy for a static magnetic 
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Figure 7: E e ff( 2 +±)/d 2 as a function of d for B m = 5 and rrif — 1 for the Gaussian magnetic 
field of Eq. (2). The discrete points correspond to our numerical results and the continuous 
line 0.448 - 0.034/d 2 to derivative expansion. 



field is given by 

^e//(2+l) = £(2+1) + ^(2+1) ( 42 ) 

where 

3 

£f ( 2°!n = / d 2 x^ / + °° 4(coth(s) - -) e- sm2 f /B ds (43) 

1 ' J 87T2 Jo S2 S 



For the special case of static magnetic fields, the above formulas were obtained in references 

In Fig. ^ we have plotted the effective energy E e jf( 2+ i)/d 2 as a function of d, for fixed 
magnetic field strength B m = 5 and rrif = 1 for the Gaussian magnetic field B\. We see, 
as expected, that in the homogeneous limit (d ^> l/y/B m ), our results agree quite well 
with those of the derivative expansion. In the small d limit (d \j\J B m ) , the derivative 
expansion fails to approximate the effective energy. See the deviations in Fig. ^|for d < 0.8. 

The case of 3 + 1 dimensions has not been included, as it is similar with the case of 2 + 1 
dimensions. 
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7 The strong magnetic field (B m ^> ml) 

In order to study the strong magnetic field (B m ^> rrij) behavior of the effective energy, we 

set A(fe) = G(kB m 2 , B m d 2 ) and make the change of variable y = kB m 1 / 2 in the relations 
(26) and (24) 

1 

/-? ^ /' - | - 00 ni 

E "'^=^L ^wJ GiM, - c)dv (45) 

where B m is a characteristic magnetic field strength defined as B m = Q/ird 2 and $ is the 
magnetic flux. 

For B m ^> m 2 we find the asymptotic expressions 

i 

^e//( 2+ i) = ^j o {G(y, B m d 2 ) - c) dy (47) 

B / \ /•+°° / \ 

£ e //(3+i) = ~ In (B ro /mj) I y (G(y, B m d 2 ) - c) dy (48) 

In the case of 3 + 1 dimensions, if we take into account the equation (13), we obtain an 
asymptotic formula for the case of strong magnetic field. 

£e//(3+i) = In (Bjmfj j o d 2 xB 2 (x) (49) 

The above equation is more general than Eq. (34), which was obtained in the limit rrifd — > 
for fixed 4>. 

For fixed rrifd, Eq. (49) yields an asymptotic formula in the case of large magnetic flux 

I r+oo 

£e//(3+i) = "24^2 ln ^ J d 2 xB 2 (x) (50) 

From Eq. (50) we obtain E e ff( 3+1 ^ = g^-0 2 ln(20) for the magnetic field of Eq. (3) and 
Beff (3+i) = 12 ~^ 2 4> 2 ln(20) for the magnetic field of Eq. (2). We see that for large the 3 + 1 
dimensional effective energy is proportional to —(f) 2 ln <p. Also, this analytical result explains 
the weak dependence on (f> of the £^ e //( 3+ i)/0 2 for the magnetic field of Eq. (3), that was 
observed numerically in reference H. 

In the case of 2 + 1 dimensions, the asymptotic formula for <p — ► +oo is obtained from Eq. 
(45) by setting mf = 0. In Fig. we have plotted the effective energy as a function of <fi for 
d = 1 and mj = 0, for the magnetic fields of Eqs. (2) and (3). Also, in the same figure we 
see that the curves wo<f) 3 ^ 2 — w\4> and vo(f) 3 / 2 — v 10 1 / 2 , which correspond to the magnetic fields 
of Eq. (2) and (3), fit our data very well. Thus, for large the 2 + 1 dimensional effective 
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energy is proportional to 3//2 . Also, we observe that the coefficients w = 0.413 ± 0.008 and 
v = 0.2769 ± 0.0005 agree closely with the coefficients obtained from the formula 



E 



e//(2+l) - 



C(-i/2) 



(51) 



which are 0.416 for the magnetic field of Eq. (3) and 0.2772 for the magnetic field of Eq. 
(2). This formula is obtained from the leading term of the derivative expansion (given by 
the Eq. (43)) if we set m/ = 0. Our numerical results show that the formula (51) gives the 
asymptotic behavior for large 4>, in the case of 2 + 1 dimensions. Note that the formula (51) 
does not include all the cases of strong magnetic field B m = 2<p/d 2 . In the case of small rrifd 
and fixed <p it is not correct. In this case the effective energy is proportional to 1/d as shown 
by Eq. (27) in section 4. 



1 - 



'magnetic field of Eq. (3)' 
0.413 phi A (3/2)-0.1 1 2 phi 
'magnetic field of Eq. (2)' 
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Figure 8: E e ff(2+u as a function of for d — 1 and rrif = for the magnetic fields of Eqs. 
(2) and (3). The discrete points correspond to our numerical results and the continuous 
lines to the best fit curves of the form w (/) 3 ^ 2 — Wi<f> for the magnetic field of Eq. (3) and 
Vo4> 3 / 2 — Vicj) 1 / 2 for the magnetic field of Eq. (2). Our numerical results for the coefficients 
are rounded to three decimal points. 



8 Conclusions 



We presented a numerical study of the fermion-induced 3 + 1 and 2 + 1 QED effective action 
in the presence of any static cylindrically symmetric magnetic field with finite magnetic 
flux. We used a simple novel approach for the fast computation of the effective action. The 
numerical work was greatly facilitated by the observation that the function A(k) — c (see Eq. 



17 



(20)) tends rapidly to zero for large k. The integrals (27) and (28) for the effective energy 
for 2 + 1 and 3 + 1 dimensional QED respectively are rendered convergent due to this fact. 

The function A(fc) — c is the only quantity we need to know in order to compute the 
effective energy (see Eqs. (27) and (28)). We have performed computations for eight mag- 
netic field configurations (see section 1). The form f\ of the function A(k) — c is similar for 
all these fields. This means that our conclusions for the effective energy are quite general 
and are not valid just for the magnetic fields of Eqs. (2) and (3). 

For e 2 /4ir = 1/137, our numerical work shows that for the families of magnetic fields we 
examined, there is no radius d of the magnetic flux tube that minimizes the total energy for 
fixed magnetic flux. The main reason for this is the small contribution (of the order of one 
per cent) of the effective energy to the total energy. An exception, where the contribution is 
not small, is the case of 3 + 1 dimensional effective energy for small d, as may be seen from 
Eq. (36). 

In section 5, we observed that our numerical results for E e ff/d 2 , in the case of the 
magnetic field of Eq. (3), converge slowly to an asymptotic value given by the Schwinger 
formula for a homogeneous magnetic field. In addition, we showed numerically that this 
asymptotic behavior for large d is of the form ao — a\/d. The case of the Gaussian magnetic 
field of Eq. (2) is different, as the convergence to the asymptotic value is quite faster than 
the case of the magnetic field of Eq. (3). The asymptotic behavior of E e ff/d 2 ~ b Q — bi/d 2 
in this case, as we see in Fig. [F], is the same with that given by the derivative expansion. 
Generalizing, we expect that this behavior would be valid for all smooth magnetic fields of 
the form we examine. In the case of discontinuous magnetic fields, for which the derivative 
expansion fails, we expect a different asymptotic behavior of the form: jd. 

We studied the case of large magnetic flux <fi (strong magnetic field) when the spatial 
size of the magnetic flux tube is fixed. We showed that the effective energy for large <fi is 
proportional to — (ft 2 In <p in the case of 3 + 1 dimensions and to 3 / 2 in the case of 2 + 1 
dimensions. Especially in the case of 3 + 1 dimensions, we derived an explicit formula for 
large magnetic flux (0 ^> 1). 

It would be interesting, if for large <j) the effective energy in the presence of an inho- 
mogeneous magnetic field would become greater than the classical energy. In view of the 
above-mentioned results we should exclude this possibility at least for the class of inhomo- 
mogeneous magnetic fields of the form of a flux tube. An exception may be the case of 3 + 1 
dimensions for very large <fi. However, for very large the asymptotic behavior — <p 2 In <p is not 
reliable, as the contributions of the higher loop corrections (two-loop corrections, three-loop 
corrections etc) to the QED effective energy cannot be assumed small. 
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A APPENDIX: Effective energy for massless fermions 

In this appendix we derive an expression for the effective energy for massless fermions, since 
Eq. (14) contains a singularity at mj = 0. It is remarkable, that this singularity does not 
exist in the unrenormalized effective energy of Eq. (12). It arises due to the renormalization 
procedure used to remove the ultraviolet divergent part of the effective energy. However, 
the above-mentioned singularity does not appear, if instead of the renormalization condition 
11(0) = we use the off-shell renormalization condition Il(— M 2 ) = 0. 

Setting rrif = in the unrenormalized effective energy of Eq. (12) we find 

£e//(3+D = 

M M 

The ultraviolet divergent part in the above equation is removed if we include the following 
counterterm in the QED Lagrangian. 

Scount = (Z 3 -I) J d A x{-- A F, v F^) (52) 

where 

Z, = l-—A-H^-) + -\ (53) 

This counterterm is defined uniquely by the off-shell renormalization condition Il(— M 2 ) = 0. 
In this case the renormalized effective action for massless fermions is 

M M 
In a similar way with that of section 3, we obtain 

E { ;; f % +1) = j 2 ( Cl (0) HdM) + cM) (55) 

where 

ci(0) = ^ Jf y(F(y, 0) - c)dy = ^ J B\x)d 2 x (56) 

c 2 (0) = -^ /' +OO yln(y)(F(y,0)- C )dy-|c 1 (0) (57) 

7T Z JO 6 

The total energy is given by the equation 

E to t = E dass + E e ff 

l r l 



?p 2 



J B\x)d 2 x + — (ci(0) \n(dM) + c 2 (0)) (5? 



where tu is the charge of the fermion defined at the scale M. The transformation law for 
the charge cm is given by the equation 

1 1 1 i ( M '\ 

ln XT (59) 



e 2 M , e 2 M 6tt 2 \M / 

which also guarantees that the total energy is independent from the renormalization scale 
M. 
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B APPENDIX: The appearance of metastable states 
for large magnetic flux </> 




Figure 9: 8^ a {k) as a function of k for I = 1, a = 1, d = 1 and = 50 for the Gaussian 
magnetic field of Eq. (2) 

In this appendix we examine the appearance of metastable states for large magnetic flux 
cf). Note that in the large magnetic flux <fi = B m d 2 /2 limit we distinguish two cases: a) a 
large magnetic field strength B m keeping d constant and b) a long range d of the magnetic 
field keeping B m constant. 

In Fig. |9], we see that the function Si yS (k) exhibits jumps at a finite number of values 
of k. Every jump increases the phase shift by 7r and corresponds to a metastable state of 
the electron. It is remarkable that these jumps occur in a small but nonzero interval Ak 
at which the density of states appears as a very sharp peak. For the calculation of the 
phase shifts in this case we have used the differential equation (75) of appendix B, which is 
appropriate for large (p. For the plot we calculated the phase shifts with a step equal to 0.1. 
It may be noted, incidentally that our numerical results for the phase shifts at k = are 
in agreement with the Levinson's theorem [20|. Levinson's theorem, for a two dimensional 



Schrodinger equation and for a potential with asymptotic behavior 0(l/r 2 ) for large r, has 
been investigated in reference ||21||. Thus for our special case 



r f 7T + £ (\l\ — \l — (j>\) if a zero mode exists ,„„ x 

^ (0) = \ f(|/|-|/-0|) otherwise (60) 

The metastable states are due to the potential well form of the effective potential fz, CT (r) + 
(/ 2 — l/4)/r 2 as shown in Fig. jT^. It is interesting to compare with the WKB approximation 



to the eigenvalues k n of the metastable states as roots of the equation 




rh [ \ l 2 \ 1 

dr=(n+- 



k*-v liS (r)- — \ dr = (n + -)7r n = 0,l,2... (61) 
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Figure 10: vi t(T (r) + - as a function of r for I — 1, a — 1, d — 1 and = 50 for the 
Gaussian magnetic field of Eq. (2) 

where a and b are the turning points. In Table [I] we see that our numerical results agree 
closely with those of WKB approximation. Note, that the electron metastable state energies 
corresponding to the eigenvalues k 2 n , are Jk 2 + m 2 f . 



n 


WKB approximation k n 


Numerical computation 


1 


13.721 


13.7-13.8 


2 


19.051 


19.0-19.1 


3 


22.881 


22.8-22.9 


4 


25.854 


25.8-25.9 


5 


28.200 


28.1-28.2 


6 


29.997 


29.9-30.0 



Table 1: We compare the energies of the metastable states as they are obtained from the 
Fig. |8| with the WKB approximation results for d = 1, I = 1, a = 1 and = 50. 

C APPENDIX: The role of the Landau levels for the 
effective action 

Taking into account the Eqs. (9) and (14), it is obvious, that the effective energy is deter- 
mined uniquely by the spectrum of the eigenvalue equation (7 m D m ) 2 \l/{ n }(x) = E^ n ^^(x). 
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In the section 5, it was shown numerically, that the effective energy of a magnetic flux tube 
of the form of the Eq.(3) tends (for d ^> 1/ y/B m ) to the result calculated from the Schwinger 
formula for a homogeneous magnetic field. One may think, that this is unexpected, as the 
spectrum of a magnetic flux tube is continuous (with zero modes if > 1) whereas the 
spectrum of a homogeneous magnetic field, has a different structure- it is discrete and con- 
sists of the well known Landau levels. According to the appendix C these two spectra are 
related in the limit of large d. Indeed, if in the analysis of the appendix C, instead of the 
Gaussian magnetic field of the Eq. (2), we have used the magnetic field of Eq. (3) we would 
find metastable states with energies almost identical with those of Landau levels (we have 
checked this numerically). Thus, the continuous spectrum of the magnetic field of Eq. (3) 
approaches, for d ^> 1/ y/B m , to the discrete spectrum of a homogeneous magnetic field, in a 
straightforward way: the each metastable energy level approaches a correspondingly Landau 
level and at the same time more metastable states are being added at the high energy end. 



D APPENDIX: Numerical calculation of the phase shifts 

In this appendix we describe two ways for the numerical calculation of the phase shifts. 

We will adapt the method which was presented by E. Farhy et al in f22|, ^| to our 
problem. We consider two linearly independent solutions, u^ ls (r), s {r) of the radial 
equation 

d 2 I 2 - - 

dr z r z 

with asymptotic behavior e ±tkr as r tends to infinity. 
We can put these solutions into the form 

4l(r) = e^My/fHPikr) (63) 
v$ t ,(r) = e^'^y/rH^ikr) (64) 

where H^ix) and h\ 2 \x) are the Hankel functions of the first and second kind respectively. 
The complex function /3/ jS (/c, r) tends to zero as r tends to infinity. 

The scattering solutions «fc,/ )S (r) of the radial equation (62) satisfy the following boundary 
conditions: 

u*,i,-(0) = (65) 

7T 111 

u k ,i, s ( r — > +°°) ~ cos(£t - - - — + 5i yS (k)) (66) 
The solutions Ukj, s { r ) are expressed as a linear combination of the solutions (63) and (64): 
u k ,i,s{r) = e iSl >°e il3l "y/rHl l \kr) + e- i5l <°e~ ip t°^Hj; 2) (kr) (67) 
Imposing the boundary condition (65) to the solution (67) we obtain 

SUk) = -Re(3 l>s (k,0) (68) 
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Substituting the solution (63) into the equation (62) we obtain the differential equation: 

ifl^k, r) + 2ik qi {kr)P' ljS (k, r) - {{3' l>s {k, r)) 2 - v l)t (r) = (69) 

where 

qi(x) = ^{HViH{%))} (70) 
The functions (3i )S (k,r) should satisfy the following conditions: 

A,.(*,+oa) = (71) 
#,„(*, +oo) = (72) 

Because the potential tends to zero very slowly 0(l/r 2 ), we begin the numerical integration 
of the differential equation (69) from a large number r max towards zero, with the conditions: 

Pi, s {k,r max ) = (-2/0 + 2 ) — (73) 

2k r r: 



max 



ft, s (k,r max ) = -(-2/0 + 2 )-—- (74) 

ZK r max 

The above asymptotic behavior of the functions f3i s (k, r) can be obtained from the equation 
(69). 

Another way to compute the phase shifts is by solving the differential equation 

— — = ~2 r Vl ^ r ^ Jl ^ kr "> cos diA k > r ) ~ N i(kr) sin5 /jS (A;,r)] 2 (75) 

with the boundary condition Si }S (k, 0) = where Ji(x) and Ni(x) are the Bessel and Neumann 
functions respectively. The phase shift is given by the limit 

5 ltS (k) = lim 5 ljS (k,r) (76) 



This method was formulated by F. Calogero in reference |24| for the 3-dimensional case. 
It was obtained for two dimensions in reference f2"5| . 

For the numerical computations of the phase shifts we have used mainly the differential 
equation (69) because it gives results faster and more accurate than the equation (75). 
However, this differential equation is stiff for large values of and small values of I. In that 
region of and / we have used the Eq. (75). 

We also compared our numerical results with those of WKB approximation (22) and we 
found good agreement for large I or large k. 
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